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<N ■ Abstract 

We develop a systematic approach to determine the large \p\ behavior of the momentum-space 

Oh: 

wavefunction, (p(p), of a one-dimensional quantum system for which the position-space wavefunc- 
tion, ip( x )> has a discontinuous derivative at any order. We find that if the kth derivative of the 
potential energy function for the system has a discontinuity there is a corresponding disconti- 
nuity in il)( k+2 \x) at the same point. This discontinuity leads directly to a power-law tail in the 

> 

t^J- ' momentum-space wavefunction proportional to l/p k+3 . A number of familiar pedagogical examples 

^1" ' 

CN ■ are examined in this context, leading to a general derivation of the result. 
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I. INTRODUCTION 



Much of the content of traditional courses in quantum mechanics consists of solving 
the time-independent Schrodinger equation in position space and applying the appropriate 
boundary conditions to find the energy eigenstates, ip n (x). 

Student understanding of the connections between the potential energy function of the 
system and the behavior of ip n (x) is becoming an increasingly important part of the un- 
dergraduate curriculum. In quantum mechanics, as in architecture and industrial design, 
it is true that ". . . form follows function . . . "— and the detailed behavior of ip n (x) (both its 
magnitude and local variation) are strongly correlated with the behavior of the potential 
energy, V(x). Analyses focusing on these connections can appear as early as modern physics 
courses (at the level of sketching wavefunctions) up through formal implementations of the 
idea via approaches such as the WKB method. Many textbooks and articles^ - - compare the 
quantum mechanical probability densities in both position- and momentum-space, |?/> n (a;)| 2 
and \4> n {p)\ 2 , and their classical analogs. 

Although some of the standard textbook level examples, such as the harmonic oscillator, 
give solutions which are infinitely differentiable, several of the most familiar model problems 
are based on potential energy functions that have discontinuities in some derivative, are 
discontinuous themselves (such as the step potential), or are singular (the 5-function and 
infinite well.) 

All one-dimensional potentials must give solutions for which i/j n (x) is continuous and for 
which at least the expectation value of p 2 (necessary for the Schrodinger equation) is well- 
defined. Higher-order derivatives of ip n {x) can be discontinuous (or singular), implying that 
expectation values of higher powers of p are not defined and repeated application of the 
differential momentum operator, p, will "find" such discontinuities in position space. 

In momentum space the expectation value of powers of momentum is given by 



and questions related to the evaluation of such average values will necessarily be tied to the 
large \p\ behavior of </> n (p) and whether it leads to convergent integrals. A natural question 
is to what extent the continuity properties of V(x) are reflected, directly or indirectly, in 
the asymptotic behavior of <fi(p), which is the subject of this paper. In other words, can we 
tell from the form of V(x) how <p n {p) behaves for large values of \p\l 




(1) 
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We will demonstrate, first through examples, and then by a more formal derivation, that if 
V(x) has a discontinuity in its kth derivative, V^(x) at x — a, then there is a corresponding 
discontinuity in ip( k+2 \x), with 

^+2) (a +) _ ^(fc+2) (a -N = ^ ry(fc)( a +) _ yW( a -)l ^( a ). ( 2 ) 

at 

The leading term in the large |p| expansion of <p(p) is directly related to this generalized 
"kink" and is given by 

m = H) « e -.w,i«^£^ A) M , (3) 

so that expectation values through (p 2 ( fc + 2 )) are well-defined. The use of Eq. (j2]) in Eq. ([3]) 
demonstrates how the form of cf>(p) depends on V^ k \ If the wavefunction vanishes at the 
cusp, if) (a) = 0, then the leading-order term will be 0(h/p) k+4 . In this language, singular 
potentials (such as the 5-function or any function involving an infinite wall/barrier) corre- 
spond to k = — 1, and discontinuous potentials such as the step potential have k = 0. The 
symmetric linear potential, defined by V(x) oc \x\ (with a cusp at x — 0), which is discussed 
in Sec. IVHt corresponds to k = 1. 

We will follow a pedagogical approach, in many ways following the path we took in 
exploring the question. We first examine the connection between the singular nature of 
V(x) and the large \p\ behavior of <p(p) for two familiar textbook-level systems, the 5- 
function potential (in Sec. [TT]) and the infinite well (in Sec. IIIII ) In Sec. |TV]we discuss the 
same issues for a single infinite wall potential, the quantum bouncer. We then illustrate 
in Sec. |V] the simple intuitive example that first led us to the systematic expansion of 
4>(p), leading to a formal general solution in Sec. I VII We include a final exemplary case in 
Sec. IVHl which illustrates a special circumstance in which the general result requires more 
careful interpretation. We then review our results, conclusions, and suggest further avenues 
of study. 

II. SINGLE (5-FUNCTION 

A simple model system with a singular potential is the single attractive 5-function po- 
tential, defined as 

Vg(x; a) = — g5(x — a). (4) 



We integrate the Schrodinger equation over the interval (a , a + ) = (a — e, a + e) and find 
the discontinuity condition on the wavefunction 

^(a + )-^(«-) = -^(a), (5) 

which can be used to determine the energy eigenvalue condition. There is a single bound 
state solution, given by 

V„(aO = V^oe-^l'-l, (6) 
where i^o = mg/h 2 , with bound state energy 

The expectation value of the potential energy is given by 

(HV(xMo) = -g\Ma)\ 2 = -gK = = -2 , (8) 

which implies that the average value of the kinetic energy is 

(MT^o) = E - (MV{x)\^) = (9) 

Equation (jUJ) can be confirmed by using either of two equivalent expressions for the expec- 
tation value of the kinetic energy 
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dx 



dx, (10) 



where an integration by parts was used to obtain the second equality from the first. From the 
second form in Eq. (flOl) . which requires only the first derivative of ip(x), we have |t/>o(o;)| = 
Ko4>o (x), giving 

(^oiT|^o) = (11) 
The first form in Eq. (Tl0|) can also be used if we note that 

d 2 ipo{x) 



K& Q (x) - 2K ^K 5(x - a), (12) 



dx 2 

where the important ^-function contribution arises from differentiating the discontinuous 
ip'(x) at x = a, using the relation Q'(x) = 5(x) for the Heaviside function. This result is 
significant because it confirms that further derivatives of ip(x) are not well defined, so that 
expectation values of higher powers of p are not calculable. 



The corresponding momentum-space wavefunction is given by 

* w ^£> w M = )/? fe) e ~ w * (13) 

where po = fo-Ko = mg/h. Standard integrals then give 



roc 



<V>ob1V>o)= / p1</>o(p)rdp = Po, (14) 

J — oo 

giving (T) = +(hK ) 2 /2m as expected. 

The large \p\ behavior of the momentum-space wavefunction is given by \<p(p)\ oc 1/p 2 
which implies that expectation values of powers of p higher than 2 will not lead to convergent 
integrals. This behavior is the first hint of the connections between the continuity behavior 
of if)(x) and the large \p\ behavior of 4>{p). 

As an aside, we note that we need to calculate expectation values only of even powers of 
p, because for stationary state solutions of bound state systems, the energy eigenfunctions 
can be put into purely real form, which implies that \<f>(— p)\ 2 = |0(+p)| 2 and expectation 
values of odd powers vanish. The most obvious example is that (i) n \p\i/) n ) = (cj) n \p\(j) n ) = 
for any such state, corresponding to the fact that the particle is equally likely to be found 
moving to the right (p = +\p\) or the left (p = — |p|). 

An alternative derivation of <p(p) which also gives the form in Eq. ffl3|) involves Fourier 
transforming the Schrodinger equation directly into momentum spaced by multiplying by 
exp(—ipx/h) and integrating over all space, which gives 

f^ (p) - -£= / + °° S(x - a) Mx) e~ tpx/h dx = -\E q \Mp), (15) 

or 

Up) = — ^ e~^' h (16) 

This form is useful because it allows an easy generalization if multiple 5-function potentials 
are present. If 

N 

V(x) = -^2g,i5(x - di), (17) 

8=1 

we immediately have that 

n N 



<j>{p) 



2m 1 
^h(p 2 + 2m\E\) 



Y,^^y~ ivailh - (18) 



i=i 
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The presence of multiple singularities in the potential (and multiple cusps in the wavefunc- 
tion) allows for interference between the exp(—ipai/h) phases, but still gives the overall 1/p 2 
behavior for large \p\. 

Looking forward to our general result, we examine the \p\ — > oo behavior of Eq. ( 1T5|) in 
a slightly different way. If we use the connection in Eq. (J5J), we find that the large \p\ limit 
can be written in the form 



This result is important because it shows that the large p momentum-space wavefunction 
can be written in a way that depends only on the properties of the cusp in the wavefunction 

clt OC —— (X. 

Although derived for this case, we will find that the result in Eqn. ( Tl9l) is the first non- 
trivial term is a systematic expansion of the leading-order large \p\ behavior of <f)(p) for which 
the corresponding ip(x) has a discontinuity in some derivative. 

III. INFINITE WELL 

For an infinite well potential defined in the region (0, L), the normalized wavefunctions 

are 



with n = 1,2,3,.... The spatial derivatives necessary to evaluate the kinetic energy using 
either form in Eq. f fl0|) are given by 




(19) 




(20) 





(21) 



and 



dx 2 




y| s . n ^ + ™^ [s{x) + { _ 1)n+l5{x _ L)]t (22) 



Using either form of Eq. (flOj) , we find that 



(if> n \f\i> n ) = —(^ n \p 2 \^n) 



H 2 n 2 TC 2 
2mL 2 




(23) 
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where p n = nnfijL. We also note that expectation values of higher order powers of p are 
ill-defined. 

The momentum-space wavefunction for the infinite well is given by 




The expectation value of p 2 can be evaluated via 

/+oo 
P 2 \MP)\ 2 dp = P 2 n (25) 
-oo 

using standard integrals, symbolic manipulation software, or by contour integration tech- 
niques. Here again, higher powers of p are not well-defined expectation values because 
\<f)(p)\ ~ l/p 2 - We also see oscillating behavior as expected from the interference of contri- 
butions from two singular potentials (the two walls). 

The large \p\ behavior of Eq. ( I24bl) is consistent with the simple form in Eq. ( fi~9l) . using 
the contribution from x = and x = L and the values 



V4(0")=<(L+) = (26) 




The form in Eq. (fl9|) has been found in the context of two simple systems, ones with 
singular potentials added to otherwise free particle systems. We next turn to a nontrivial 
system to which an impenetrable boundary has been added, to further explore the generality 
of the result in Eq. f fT9|) . 



IV. QUANTUM BOUNCER 

Another pedagogically familiar system which includes a single infinite wall barrier is the 
quantum bouncer,- defined by the potential 

Fz (z > 0) 

V(--) = { • (29) 

oo (z < 0) 



This potential has received renewed interest as the simplest model for recent experiments 
showing evidence for quantized energy states of neutrons in the Earth's gravitational field.— 



• 8-1(1 



Other recent applications of this model to physical problems are discussed in Refs. 

The properties of the Airy function solutions have been re-examined in the light of this 
renewed interest with new analytic results,—"— using identities which appeared some time 
ago.— For example, the normalized wavefunctions are, 

jN n Ai(y-Q (y>0) 
ip n (z) = < (30) 

(y<0) 



where 



N - = 7m=ur (31) 



The relevant dimensionless quantities are 

z ( h 2 \ 1/3 

y-- P and '=UsfJ • (32) 

The — ( n are the zeros of the well-behaved Airy function, Ai(z), and the energies are given 
in terms of them as 

h 2 F 2\ 1/3 



E n = £o(n where £ = Fp= — — . (33) 

\ 2m J 

The corresponding classical position-space probability density is 

P n CL) {x) = 1 (0 < z < An) (34) 

2 v /A n (A n - z) 

and zero elsewhere. The upper classical turning point is given by ( n £o = E n = FA n , or 
A n = p( n - The classical momentum-space probability density is given by 

1 



-Qn<P< +Qn) 



P^\ P ) = { 2 Qn ' (35) 

otherwise 



where 



( n £ = E n = & or g„ = - v / Cn- (36) 
2m p 

We note that the classical and quantum probability densities for both position and momen- 
tum for a closely related system, the symmetric linear potential 2 (discussed in Sec. IVIIj) . 
have been compared. 
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The momentum-space wavef unctions can be obtained numerically by using the ip n {z) in 
Eq. (130]) and the definition of the Fourier transform, 

i r+oo 

MP) = -f= / Mz) e- ivzth dz (37a) 
V^nti Jo 

r+oc |- ( pz\ {PZ\~\ 

= v^nJo Mz) \- cos th) ~ lsin Vh)\ dz (37b) 

= ^{p)-i^(p). (37c) 

We plot the contributions of \4>n (p)\ 2 , \4>n (p)\ 2 , and their sum, each compared to 
Pi CL \p), in Figs.^a), (b), and (c). 

Note that the contribution of the imaginary component of <p{p) shown in Fig. 1(a) is sim- 
ilar to many of the visualizations of comparisons between classical and quantum mechanical 
probability distributions, where the quantum result oscillates about the classical prediction, 
consistent with WKB-type approximations. In this case, the imaginary component provides 
half of the total probability, in that "locally averaged" sense. The real part [see Fig. 1(b)] 
has similar oscillatory behavior, but with a slightly different structure. The combination of 
the two gives a much smoother approach to the classical "flat" momentum distribution. 

More importantly, we see that the real component of the Fourier transform (the one 
affected most directly by the infinite wall) extends much further into the classical disallowed 
region of momentum space, hinting at the expected power-law 'tail'. We can fit the large 
\p\ tails of 0n RE \p) for various values of n and find the simple result (p and H set equal to 
unity) 

P 2 |0i RE) (p)I^O.4 (p-Kx>), (38) 

independent of n. To compare this result to the prediction of Eq. (1T9T) . we note that ^(O - ) = 
0, and from the normalized wavefunction in Eq. (1301) . we find that 

which is independent of n, so that in dimensionless form 

P 2 \4>1 RE) (P)\ ~-L^ 0.3989, (40) 

V Z7T 



which is another important confirmation of the simple result in Eqn. ( fl9l) for singular poten- 
tials. (We note that a similar numerical evaluation and subsequent fitting of $n\p) finds 
that it scales as 1/p 5 . We will discuss this result in Sec. IVIII ) 



V. HYBRID EXAMPLE: HINTS OF THE GENERAL SOLUTION 



We next consider potentials that are discontinuous, but not singular. The most familiar 
example is a step potential, defined by 

V s (x;a) = V Q(x -a) (41) 

where Q(x) is the Heaviside function. For a finite well (FW) we have can write 

(x < a) 

V FW {x; a, b) = -V s {x; a) + V s {x; b) = l-V Q (a < x < b) ■ (42) 

(x>b) 

The boundary conditions for such discontinuous potentials are that both if) and if>' are 
continuous across such a step,— and thus we might expect a qualitatively different behavior 
of 4>(p), due to the poor behavior of if)" at a step boundary. 

To pursue this question and to allow for a more systematic study of the large \p\ behavior 
of 4>(p), we consider the hybrid case of an attractive 5-function potential combined with a 
single step potential. For definiteness, we consider 



V(x) = V s (x;0) + V s (x;a), 



(43) 



which is a slight generalization of a single 5-function potential interacting with an infinite 
wall, as discussed by Aslangul.— We note that for sufficiently large positive V , the single 
bound state is no longer supported, and if V < 0, the possibility of tunneling can also 
preclude a stable bound state. Although the study of what range of Vq values support a 
bound state is an interesting question in itself, we assume that Vq is such that there is one, 
with energy E = — \E\, and focus on the behavior of the corresponding <f>(p). We do not 
need to find a normalized solution in detail because we will focus only on the nature of the 
wavefunction (and any discontinuities) and the locations of the singularity and discontinuity 
in the potential. 

The (un-normalized) solutions for each region can be written as 



if>(x) 



Ae +Kx (x < 0) 

Be~ Kx + Ce +Kx (0 < x < a) 
De~ Qx (x > a) 



(44) 
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where 



2 2m\E\ 2 2m(\E\+V ) 

K = — — — and Q = — 

nr n 2 



(45) 



The boundary conditions which are imposed by the singular 5-function and discontinuous 
step potential are 



ip(x = 0) 
i/)(x = a) 
ij)'(x = a) 



A = B + C 

Be- Ka + C +Ka = De- Qa 
- KBe~ Ka + KCe +Ka = -QDe~ Qa . 



(46a) 
(46b) 
(46c) 



The corresponding Fourier transform is 

1 f A 



<f>(p) 



+ 



B 



C 



V2nh [K -ip/H K + ip/h K-ip/h_ 
+ J— { Be ~ Ka + Ce+Ka + De ~ Qa \ e -W 



y/2nfi \ K + ip/H + K - ip/h + Q + ip/H } 
<t>s{p) + <t>s{p). 



(47a) 
(47b) 



where we have separated the terms related to the singularity and discontinuity. 

Because we are interested in the large \p\ behavior, we can systematically expand each 
term in inverse powers of p. For example, we find that for the S-f unction contribution 

1 



f(p) 



-i(A -B-C)- + (AK + BK- CK)% + 
p p 2 



(48) 



The first term vanishes because of the continuity of ip at the origin, and the second one 
is consistent with the general form for singular potentials in Eq. ( [19]) (with the singularity 
located at x = so that exp(—ipa/h) — 1.) 

For the term arising from the single-step potential, we have 



V2nh 



i(Be~ Ka + Ce 



-Ka 



De 



-Qa 



h 



h 2 



- + (-Bke~ Ka + CKe +Ka + Dqe~ Qa ) — 



+ i(-BK 2 e 



2-Ka 



h 3 



CK e + + DQ e~ Qa )— + 



-ipa/h 



(49) 



where we have carried the expansion to one higher order. The first two terms vanish because 
of the continuity of ip and ip' given by Eqs. (I46al) and (146bj) respectively. The remaining 
non-vanishing higher order term can be expressed as 



V2vr^ P A 



(50) 
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Note that for discontinuous potentials, we can derive a relation similar to Eq. (jSJ) for a 
step-potential of the form in Eq. (jHJ). By differentiating the Schrodinger equation once, and 
then integrating over the range (a - , a + ) we find that 

r{a + )-r{a-) = -^4>{a) (51) 

because V s '(x; a) = V 6(x — a), which can be used to simplify Eq. fl50|) . 

Motivated by these results, we have confirmed that the expression in Eq. ( )50|) gives the 
leading large \p\ behavior of <f>(p) for a variety of model systems involving discontinuous step 
potentials, including the finite well of Eq. ( H2|) . the asymmetric finite well (with different 
barrier heights on either side), and other combinations. 

More importantly, if we use this problem as a template, we find that we can express the 
series expansion for each term, 4>s,s{p), in Eqs. (HHj) and (f49|) . in the common form 

^(p) = -^=X) r *(p)« ( 52 ) 



where 



2nh k=0 



T^p) = -i \4){a + ) - if)( a -)] e~ ipa/h - (53a) 

p 

75 (p) = - W{a + ) - ^\a-)] e- ipa / h ^ (53b) 

%(p) = +i ty"(a + ) - ^"(a-)] e"^/' 1 ^ (53c) 

7i(p) = + [r{a + ) - r'(a~)] e~^ h ^ (53d) 



%{p) = (-2)™[^ (n - 1) (a+) - ^ n - l \a-)} e~ ipa/h ( - ) . (53e) 



In each case, the vanishing of the leading 0(p~ l ) term is a consequence of ip(x) being 
everywhere continuous, which is required for a consistent probability interpretation. This 
condition then implies that the lowest-order term possible for <fi(p) is of order 0(p~ 2 ) for 
large \p\, which ensures that (p 2 ) always gives a convergent integral. 

VI. FORMAL SOLUTION 

To confirm the systematic expansion of 0(p) suggested by Eq. ( 153|) . we now describe an 
approach to the evaluation of the Fourier transform <p(p), focusing on the case where the 
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corresponding ifj(x) has a discontinuous derivative at some order. For simplicity, we assume 
that the generalized kink is at x = 0; the extension to any other location is trivial. 
We first write the Fourier transform as 



(f>(p) 



1 



V2ttE 
1 

1 

\phxfi 



if>(x) e~ ipx/h dx 



oo 



ip(x) COS 



[h(p)-ih(p)\ 



px 



dx — i- 



V2ttE 



ip{x) sin 



px 



dx 



(54a) 
(54b) 
(54c) 



and we will focus on the approximation of I\,2{p) for large \p\. 

Assuming that any continuity issues are localized at x = 0, we split the integral into two 
regions, adding appropriate convergence factors, exp(±er), namely. 



hip) — J ip{x) cos (Jjpj dx 



lim 
lim 



ip(x) COS 



px 



e +ex dx + 



ip(x) COS 



px 



e ex dx 



li e \p) 



(55a) 

(55b) 
(55c) 



We can first rewrite the cosine term and expand ip(x) in a series expansion (one for each 
integration region) via 



cos 



h 



-ipx/h 



) 



and 



i) {n \Q)x n 



n=0 



(56) 
(57) 



giving 



Ln=0 



V> (n) (o- 



E 

n=0 



^ {n) (0 4 
n! 



x n ^x(e+ip/h) + e x( 6 -ip/R)j rfx 



/■+00 

/ 

JO 



n \ f ,- x ( e - i P/ h ) 



(58) 



We use 0~ (0 + ) for the x < (x > 0) integrals respectively. (This type of regularization 
method is similar to that used in establishing the completeness relations of the eigenfunctions 
of the single 5-function potential.—) 



13 



Because of the convergence factors, the integrals are easily performed using 

"0 r+oo 

y n e y dy = (-l) n n\ and / y n e~ y dy = n\. (59) 



These integrals give 



+ 



1 



n=0 



n=0 
oo r . n+1 



> - ip/h) n+1 (e + ip/h) n+1 
i n+1 [l + (-l) n+1 ]fr n+1 



ty(*0(0+) + (-l)^(») (0-)] (60a) 



P 



,n+l 



[ip {n) (0 + ) + (-l)V n) ((T)] ase^O (60b) 



fi 2 fi 4 fi 6 

.^[^'(o + ) - ^'(0-)] + ^'"(o + ) - ^"(o-)] - ^ (5) (o + ) - ^ (5) ((T)] + 



(60c) 



Thus, the cosine component of the Fourier transform giving <f)(p) gives the systematic ex- 
pansion in differences of odd powers of derivatives (?// 2fe+1 )(0)) at the discontinuity, just as 
in Eqs. ( l53b|) and ( l53dl) . 

In the same manner, we can evaluate the sin(px/h) integral and find 

1 1 



-*7< e) (p) 



2i 



E 



E 

n=0 



_(e-ip/h) n+1 (e + ip/h) n+l _ 



ty,(*)(0+) _ (-l) n V (n) )(0 - )] (61a) 



n+1 



[^W(0+) + (_i)^W(o-)] ase^O 



(61b) 



= -^-[^(o + ) - ^(o-)] + 4^'(o + ) - - 4^ (4) (o + ) - ^ (4) (o-)] + • • • , 
p p 6 p° 

(61c) 

which gives the appropriate even powers of derivatives, along with the correct factors of ±i, 
reproducing the results in Eqs. ([53aj) and (J53cj). Taken together, these two results give the 
general form in Eq. (15 3 el) . 

If the discontinuity is located at another location, we can use the result that if 
i/j(x) — > ip(x — a), then the corresponding momentum-space wavefunction satisfies <fi(p) — > 
<p(p) exp(—ipa/h) with the appropriate derivatives now evaluated at a ± . 

If the potential (and resulting solutions) are well-behaved, as with the harmonic oscillator, 
then all such differences of the derivatives vanish, implying that there are no power-law tails. 
This connection is not surprising because such an asymptotic expansion will not capture 
information on more well-behaved (for example, exp(— a 2 p 2 /2)) functions. 
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As noted, if the potential has a generalized discontinuity, then repeated differentiation of 
the Schrodinger equation yields a relation between the lowest-order difference in derivatives 
and the wavefunction at the discontinuity. Examples include the relations in Eqs. §5§ and 
(|5ip . If V^ k \x) has a discontinuity at x — a, then dV^ k+1 \x)/dx oc 5(x — a) and differenti- 
ating the Schrodinger equation k + 1 times, and then integrating over the range (a — e, a + e) 
yields a relation of the form 

^(*+2)(a+) - ^ k+2 \a-) = ^ [V^(a + ) - V^(a~)] 4>(a), (62) 

which can be used to simplify the leading order term in Eq. (I53ep . In the case where ip(a) = 
(where the wavefunction vanishes at the discontinuity), the expansion will necessarily start 
with one more power of p~ x . In that case, the first non-vanishing term in the (k + l)th 
differentiation of the V(x)tp(x) term gives 

^(*+3)( a +) _ ^( fc +3)( a -) = {k + 1)^ [V { - k \a + ) - V (k \a~)] f(a), (63) 

which can be used in the % + 4(p) term in the expansion in Eq. ( 1521) . A case where this 
situation occurs is considered in Sec. IVII1 



VII. THE SYMMETRIC LINEAR POTENTIAL 



We now consider another example to confirm the higher order predictions, and to note 
one additional novel feature. A straightforward generalization of the quantum bouncer is 
the symmetric linear potential, defined by 

V(z) = F\z\. (64) 

which has been described as a "parity extended" version of the quantum bouncer— and 
shares many of the same features. For this potential, which has a cusp (discontinuity in 
V'(x)), we expect that 4>{p) oc l/p 4 as in Eq. (153dl) . at least for those states with -?/>(0) ^ 
as in Eq. ( 162]) . We note for future reference that V"(z) = 2FS(z). 

Because of the symmetric nature of this potential, the eigenstates have definite parity. 
The odd states are simply related to those of the quantum bouncer and using the same 
notation as in Sec. IIVI we have 

, ' +N n Ai(y -Cn) (V>0) 

^n\ Z ) = > (65) 



y/2 



-N n Ai(-y-( n ) (y<0) 
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with the same notation and energies as in Eqs. (1321) and (I33I) . The corresponding (properly 
normalized) even states can be written in the formal 

^+\z) = M n Ai{\y\-r ln ) where M n = \ - -. (66) 

y/2pT] n Al(-T] n ) 

The —r] n are the zeros of the derivative of the well-behaved Airy function, given by 
Ai'{-r} n ) = 0. 

For this potential we expect the leading order term in the large- \p\ expansion of <f>(p) in 
Eq. §31D to be 

4=[f(a + )-f(0]^^. (67) 
V 2,nh " V 

which depends on the third derivatives of ip. If we use the same approach leading to Eqs. (jSJ) 
and 05 ip . and more generally Eq. ( 1621) . namely integrating the Schrodinger equation (differ- 
entiated an appropriate number of times) we find a constraint on if)'", namely, 

^-ty»'(0 + )-r'(0-)] =2i^(0), (68) 



or 



[tT(o + ) - V w (o-)] = (69) 

The momentum-space solutions corresponding to even states of this potential, where ip(0) ^ 
0, are expected to scale as 

V'\^\V)\^^=^= (70) 

V27T V 2 ?7n 

when all dimensional quantities are removed. This result has an n-dependence through t] n , 
which scales as n 2//3 for large values of n. 20 (Note that the (f>n\p) receives contributions only 
from the cos(px/h) term in the Fourier transform.) 

In contrast, the odd states have a vanishing value of ip(0), and thus necessarily have an 
asymptotic \p\ dependence for (p(p) that starts at one order higher, namely proportional to 
7~ n =5 ~ p~ 5 - For this case we can also simplify the general expression by writing 

<>-)-V/'"(0+) = -^M (71) 

which is most easily obtained by repeated differentiation of the Airy differential equation 
on either side of the cusp at x = 0, as in Eq. (16"5)) . For comparison to numerically obtained 
results, we once again set dimensional quantities equal to unity, in which case the large |p| 
behavior of 4>{p) is expected to be 

?W<P)\^± = ±, (72) 
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independent of n. (In this case (fin (p) receives contributions only from the sm(px/h) term 
in the Fourier transform.) We have numerically evaluated (fin (p) for several values of n and 
present the results for representative even and odd cases (n = 11) in Fig. [2l along with the 
large \p\ approximations in Eqs. ( 1701) and (1721) and find excellent agreement. 

Note that the 1/p 5 behavior of the imaginary component of (fi n (p) for the quantum 
bouncer, which here is directly related to (fi^(p), can be understood as a special case of the 
general results in Eq. (163|) . 

VIII. CONCLUSIONS 

The ability to make connections between the potential of a system and the resulting 
wavefunctions in position space is an important part of the toolkit of any quantum mechanic. 
Being able to visualize these connections and understand them at a conceptual level is an 
increasingly important focus of pedagogy 2 - in the subject. 

The corresponding skills involved in understanding and interpreting quantum phenomena 
in momentum space are far less developed. Some simple examples in this area exist, such 
as the "flat" momentum distribution for the quantum bouncer (as in Fig. [1]) arising from 
the constant force law for that system, and some others.- 1 ^ Any new examples that provide 
tangible relations between the potential energy function and <fi(p) in a direct way are valuable. 
Our results extend the study of the momentum-space probability distribution from the semi- 
classical limit to the deeply quantum regime of large momenta, far beyond the classical 
turning points in p-space. 

The single (^-potential led us to some early intuition about the structure of the general 
result, namely the form in Eq. (fl9j) . This observation was important as it emphasized the 
likely appearance of differences of derivatives at various orders depending on the behavior 
of V(x). As importantly, in generalizing this expression to Eq. (I53el) . we noted that such 
a result would be consistent with simple dimensional analysis arguments, because for every 
higher derivative used, another power of h/p would be required to compensate. 

We suggest several possible avenues for research amenable to exploration by motivated 
undergraduates. The system in Sec. IVHI can, for example, be extended to an asymmetric 
linear potential by having different force constants for x > (F) and x < (F). The 
solutions can still be written in terms of Airy functions, but the parity symmetry is now 
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broken and all states should have <p(p) at large \p\ proportional to p~ 4 , but with coefficients 
that depend on (F — F). We can also imagine extending our work in other ways, for 
example, by looking at if the correlations found here between kinks in position-space and 
tails in momentum-space are obvious in the Wigner distribution, which is one of the canonical 
quantum mechanical formats for discussing x-p connections and correlations. 

Extensions to more realistic three-dimensional systems are also natural, such as finite 
well models of the nuclear force. The Coulomb problem for the hydrogen atom, involves a 
singular potential, and the ground state momentum-space wavefunction 

(73) 

(where p = h/a$ with a$ the Bohr radius) gives only a finite number of well-defined 
expectation values for (p k ) because of its power-law (p~ 4 ) behavior for large \p\. Such 
momentum-space solutions are of interest because 4>is{p) has been directly measured in 
scattering experiments.— 
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FIG. 1: Plots of the classical momentum-space probability distribution, Pn (p) versus p, for 
the quantum bouncer (dashed curve) from Eq. (|35p and the corresponding quantum distribution, 
I<?Kp)| 2 (solid curves) from Eq. ()37c[) . The imaginary and real parts of |</>(p)| 2 from are shown in (a) 
and (b) and the total |</>(p)| 2 in (c). The vertical dashed lines labeled =LQ n indicate the classical 
turning points in momentum-space. 
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FIG. 2: The solid curves show |0(p)| 2 for the symmetric linear potential in Sec. IVIII obtained by 
numerical Fourier transform for a typical (p)) and odd (4>n \p)) case (n = 11). The 

straight dashed lines correspond to the lowest-order non-vanishing predictions in Eqs. (|70p and 
(|72p . The vertical dot-dash line indicates the classical turning point in momentum-space at +Q n 
for positive p. Note the log-log scales. 
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